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Chapter  9 

Interference  Rejection  and  Management 


Amn  Batra,  James  R.  Zeidler,  John  G.  Proakis,  and  Laurence  B.  Milstein 


9.1  Introduction 


In  this  chapter,  we  consider  interference  suppression  in  several  different  wireless 
communication  systems.  The  first  topic  treated  is  self-interference  encountered 
among  cooperating  systems,  for  example,  the  self-interference  that  is  encountered  in 
cognitive  radio  systems  and  ultra- wideband  (UWB)  communication  systems.  Both 
single-carrier,  direct  sequence  signals  and  multicarrier  signals  are  considered,  and 
the  effects  of  the  interference  on  the  performance  of  a  direct  sequence  UWB  sys¬ 
tem  that  employs  channel  state  estimation  in  the  presence  of  the  interference  is 
evaluated. 

The  second  topic  that  is  treated  is  the  mitigation  of  narrowband  interference 
in  block-modulated  multicarrier  systems.  Two  schemes,  multicarrier  code-division 
multiple  access  (MC-CDMA)  and  orthogonal  frequency  division  multiplexing 
(OFDM),  are  investigated.  The  inherent  frequency  diversity  of  MC-CDMA  through 
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the  use  of  spreading  codes  allows  for  robustness  against  interference.  Conversely, 
OFDM  must  employ  a  signal  processing  technique  to  suppress  the  interference  due 
to  its  lack  of  frequency  diversity.  The  performance  of  OFDM  is  improved  with  the 
addition  of  forward  error  correction  (FEC)  coding  in  the  frequency  domain,  thereby 
providing  the  system  with  frequency  diversity. 

The  third  topic  is  the  suppression  of  interference  in  multiple-input,  multiple- 
output  (MIMO)  wireless  communication  systems  that  employ  multiple-transmit 
and  multiple-receive  antennas  to  increase  the  data  rate  and  achieve  signal  diver¬ 
sity  in  fading  multipath  channels.  Interference  in  MIMO  wireless  communication 
systems  usually  consists  of  intersymbol  interference  (ISI)  due  to  channel  multipath 
dispersion  and  cross-talk,  or  interchannel  interference,  due  to  the  simultaneous 
transmissions  from  the  multiple-transmit  antennas.  The  focus  of  the  section  is  on 
point-to-multipoint  (broadcast)  MIMO  systems  in  which  the  channel  characteristics 
are  assumed  to  be  known  at  the  transmitter,  so  that  interference  mitigation  can  be 
performed  at  the  transmitter. 


9.2  Self-Interference  Among  Cooperating  Systems 


This  section  is  concerned  with  the  operation  of  systems  that  have  been  specifically 
designed  to  operate  in  the  same  frequency  band  as  do  other,  preexisting,  systems. 
Such  overlay-type  operations  are  often  referred  to  as  dynamic  spectrum  allocation. 
Because  the  spectral  bands  are  the  same  for  both  systems,  some  type  of  interfer¬ 
ence  suppression  and/or  interference  avoidance  is  typically  necessary.  Further,  the 
presence  of  the  interference  makes  virtually  all  standard  system  functions,  such  as 
synchronization  and  diversity  combining,  less  accurate.  In  Section  9.2.1,  a  brief 
overview  of  waveform  design  appropriate  for  overlay  systems  is  presented,  and  in 
Section  9.2.2,  the  effects  of  overlaid  interference,  both  with  and  without  interfer¬ 
ence  suppression  at  the  receiver,  on  the  accuracy  of  channel  state  information  is 
presented. 


9.2.1  Interference  Suppression  to  Enable  Spectrum  Sharing 


A  fairly  recent  goal  is  to  allow  multiple  classes  of  users  to  share  common  spectrum 
in  order  to  increase  the  spectral  efficiency  in  whatever  band  is  under  consideration. 
This  has  led  to  technologies  such  as  cognitive  radio  and  ultra  wideband  communica¬ 
tions,  which  have  the  characteristics  of  enabling  secondary  users  to  share  spectrum 
with  primary  users,  providing  the  former  set  of  users  do  not  impose  excessive  in¬ 
terference  on  the  latter  set.  Such  scenarios  typically,  require  the  use  of  interference 
suppression  and  management  techniques. 

This  involves  the  design  and  analysis  of  spectrum  sensing  techniques  in  or¬ 
der  to  determine  where  and  when  primary  users  are  present.  It  also  involves 


9  Interference  Rejection  and  Management 


219 


choosing  appropriate  waveform  designs,  such  as  single-carrier  DS  or  multicarrier 
techniques  (e.g.,  multicarrier  DS),  and  combining  them  with  appropriate  signal  pro¬ 
cessing  schemes  to  enable  them  to  both  withstand  the  affects  of  overlaid  interference 
and  intentional  jamming  and  minimize  the  interference  they  impose  on  the  overlaid 
waveforms.  For  the  latter  case,  the  signal  processing  will  take  the  form  of  placing 
appropriate  notches  in  the  spectrum  of  the  transmitted  DS  signals  at  frequency  lo¬ 
cations  where  the  spectrum  sensing  operation  indicated  that  preexisting  narrowband 
users  were  active.  In  the  former  case,  the  signal  processing  can  again  take  the  form 
of  notch  filtering,  but  now  at  the  DS  receivers.  In  other  words,  notch  filtering  at 
the  DS  transmitters  is  designed  to  protect  the  overlaid  users  and  notch  filtering  at 
the  DS  receivers  is  designed  to  protect  the  DS  signals  themselves.  Also,  the  term 
“notching”  is  used  in  a  general  sense,  meaning  that,  for  example,  in  the  case  of  mul¬ 
ticarrier  DS,  the  notching  might  imply  simply  not  transmitting  appropriate  subcar¬ 
riers.  Note  that  the  interaction  of  channel  estimation  errors  (such  as  those  incurred 
while  performing  tasks  such  as  spectral  sensing  and  channel  state  estimation)  and 
interference/jamming  is  an  important  consideration  and  is  discussed  below. 


9.2.1. 1  Single-Carrier  Direct  Sequence 

Adaptive  MMSE  receivers:  As  is  well  known,  a  Wiener  filter  can  greatly  improve  the 
performance  of  DS  receivers  in  the  presence  of  narrowband  interference.  Based  on 
the  a  priori  knowledge  of  the  channel  and  narrowband  interference  (NBI)  statistics, 
a  Wiener  filter  is  able  to  suppress  NBI  and  obtain  a  good  estimate  of  the  channel 
realization.  For  instance,  the  frequency  response  of  a  Wiener  filter  is  presented  in 
Section  9.2.2  for  different  NBI  levels,  and  it  will  be  seen  that  the  filter  suppresses 
the  frequencies  where  the  NBI  is  located.  In  practice,  the  channel  and  narrowband 
interference  statistics  are  not  known  and  might  be  changing  with  time.  Therefore, 
the  Wiener  filter  must  be  made  adaptive.  As  shown  in  [49,  53],  the  LMS  filter  is 
very  effective  in  rejecting  NBI  when  the  ratio  of  the  NBI  bandwidth  to  the  signal 
bandwidth  is  small,  assuming  the  LMS  filter  has  had  sufficient  time  to  converge. 

Tunable  notch  filters:  There  are  many  examples  in  the  literature  that  illustrate 
the  usefulness  of  notch  filters  when  trying  to  suppress  narrowband  interference.  For 
example,  it  was  shown  in  [33]  that  the  presence  of  an  interference  rejection  filter  can 
significantly  improve  the  acquisition  system  performance  of  a  DS  CDMA  receiver. 
And  it  was  shown  in  [34]  that  in  order  to  successfully  have  a  CDMA  system  overlay 
narrowband  users,  i.e.,  to  deploy  it  in  a  manner  such  that  neither  set  of  users  caused 
excessive  interference  to  the  other  set,  it  was  desirable  to  use  tunable  notch  filters  at 
both  the  CDMA  transmitters  and  the  CDMA  receivers. 


9.2. 1.2  Multicarrier  Direct  Sequence 

In  a  multicarrier  DS  system,  multiple  narrowband  DS  waveforms,  each  at  a  distinct 
carrier  frequency,  are  combined  to  yield  a  composite  wideband  DS  signal.  Among 
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the  advantages  of  such  an  approach  is  the  ability  to  achieve  the  same  type  of  sys¬ 
tem  performance  that  a  single-carrier  DS  signal  would  provide,  such  as  diversity 
enhancement  over  a  multipath  channel;  however,  this  is  achieved  without  the  need 
for  a  contiguous  spectral  band.  It  was  shown  in  [9]  that  for  an  exponential  mul¬ 
tipath  intensity  profile  and  imperfect  channel  estimation,  multicarrier  DS  systems 
outperform  single-carrier  DS  systems. 

In  the  multicarrier  scenario,  the  receiver  first  operates  in  a  spectrum  sensing 
mode,  detecting  the  presence  and  frequency  occupancy  of  possible  narrowband 
users,  as  well  as  possible  jammers.  Then,  after  appropriate  subcarriers  are  excised, 
both  channel  estimation  and  data  detection  are  performed. 


9.2.2  Effects  of  Interference  on  Channel  State  Estimation 


One  of  the  key  needs  of  wireless  receivers  is  accurate  channel  state  information 
(CSI)  so  that  basic  functions  such  as  synchronization  and  diversity  combining  can  be 
accomplished.  In  this  section,  we  consider  the  problem  of  trying  to  achieve  accurate 
CSI  in  the  presence  of  external  interference;  we  also  look  at  the  impact  on  system 
performance  when  either  the  estimation  module  or  the  data  detection  module  at 
the  receiver  does  not  take  sufficient  measures  to  account  for  the  interference.  As  a 
specific  system  with  which  to  illustrate  the  concepts,  we  will  use  UWB. 

Consider  a  direct  sequence  UWB  system  as  described  in  [43,44].  A  block  dia¬ 
gram  of  the  low-pass  equivalent  receiver  is  shown  in  Fig.  9.1.  The  first  block  is  a 
chip-matched  filter  which  is  perfectly  synchronized  to  the  desired  waveform,  and 
where  binary  phase- shift  keying  (BPSK)  is  the  modulation  format.  The  output  of 
the  matched  filter  is  sampled  at  the  chip  rate,  and  the  samples  are  then  processed 
by  a  rake  receiver.  The  received  signal  consists  of  the  UWB  waveform,  additive 
white  Gaussian  noise  (AWGN),  and  narrowband  interference  which  is  modeled  as  a 
wide-sense  stationary  Gaussian  random  process. 


Fig.  9.1  Low-pass  equivalent  of  the  DS/CDMA  receiver.  For  the  reception  of  the  vth  data  bit, 
Wmin(0  =  (v  —  1)  A/^  +  /  and  mmax(0  =  vN  +  /  —  1.  For  the  vth  channel  estimation  period, 
Wmin(0  =  (r  —  1)  NNp  +  i  and  mmax(0  =  vNNp  +  /  —  1.  Note  that  N  is  the  processing  gain  and 
Np  is  the  number  of  pilot  symbols. 
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There  are  two  types  of  receiver  configurations  that  are  studied.  The  first  one  is 
based  upon  simplicity  of  implementation:  It  uses  a  sample  mean  (SM)  to  estimate 
the  gain  of  each  resolvable  path  used  in  the  rake  receiver  of  Fig.  9.1  and  com¬ 
bines  those  paths  that  are  resolved  using  maximal-ratio  combining  (MRC).  Note 
that  neither  the  estimator  nor  the  detector  requires  knowledge  of  the  statistics  of  the 
channel,  and  neither  of  them  make  any  explicit  attempt  to  suppress  the  interference. 
The  second  receiver  is  more  complex.  It  employs  a  minimum  mean-squared  er¬ 
ror  (MMSE)  estimator  and  a  maximum-likelihood  (ML)  detector.  Thus,  the  second 
receiver  requires  knowledge  of  the  second-order  statistics  of  the  interference,  and 
those  statistics  typically  are  unknown  at  the  receiver  and  thus  requires  an  estimate. 
In  what  follows  we  assume  that  this  latter  information  is  available  noise  free. 

A  detailed  analysis  of  this  problem  is  presented  in  [43,44].  Here,  we  just  describe 
some  of  the  results  in  the  context  of  the  effect  that  the  interference  has  in  achiev¬ 
ing  accurate  CSI.  Consider  first  Fig.  9.2.  The  leftmost  curve,  labeled  (I),  displays  the 
performance  of  the  system  when  perfect  CSI  is  available  (i.e.,  we  have  a  genie-aided 
receiver)  and  ML  detection  is  used.  This  curve  corresponds  to  interference-free  con¬ 
ditions,  and  thus  represents  a  lower  bound  on  the  performance  of  any  of  the  other 
systems  to  be  discussed.  The  three  curves  labeled  (II)  correspond  to  imperfect  CSI 
and  to  signal-to-interference  ratios  (SIR)  of  —10,  —15,  and  — 20dB.  The  receiver  is 
an  MMSE/ML,  and  we  observe  that,  because  this  receiver  has  statistical  knowledge 
of  the  interference,  it  can  configure  itself  to  suppress  that  interference.  This  will  be 


Fig.  9.2  Bit  error  rate  for  the  DS/CDMA  system:  (I)  perfect  channel  estimation  and  ML  data 
detection  in  the  absence  of  NBI,  (II)  MMSE/ML  receiver,  and  (III)  SM/MRC  receiver.  SIR  = 
—  lOdB  {solid  lines),  —15  dB  {dashed),  — 20dB  {dash-dot).  Number  of  training  bits:  50 
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illustrated  below  and  is  the  reason  the  latter  three  curves  of  Fig.  9.2  are  so  close  to 
the  interference-free  curve. 

If  we  replace  the  MMSE/ML  receiver  with  an  SM/MRC  receiver,  we  get  the 
performance  curves  shown  in  Fig.  9.2  that  are  labeled  (III).  Note  that  all  three  curves 
now  experience  an  error  floor,  and  that  the  performance  of  the  curve  with  an  SIR  of 
—20  dB  is  virtually  worthless,  notwithstanding  the  fact  that  the  system  is  operating 
with  a  processing  gain  of  127  (i.e.,  roughly  21  dB). 

All  the  estimates  used  in  the  curves  shown  in  Fig.  9.2  were  made  with  a  training 
sequence  of  50  bits.  To  see  the  sensitivity  of  the  system  to  the  number  of  train¬ 
ing  bits,  consider  the  curves  shown  in  Fig.  9.3.  As  in  Fig.  9.2,  the  leftmost  curve 
corresponds  to  perfect  CSI,  the  curves  labeled  (II)  correspond  to  the  MMSE/ML 
receiver,  and  the  curves  labeled  (III)  correspond  to  the  SM/MRC  receiver.  For  each 
of  the  latter  two  sets  of  curves,  results  are  shown  for  observation  intervals  corre¬ 
sponding  to  10,  30,  and  50  bits.  Once  again,  it  is  seen  that  the  performance  of  the 
MMSE/ML  receiver  in  the  presence  of  channel  estimation  error  is  very  close  to  the 
performance  achievable  with  perfect  CSI,  whereas  the  performance  of  the  SM/MRC 
receiver  suffers  tremendous  degradation. 


Fig.  9.3  Bit  error  rate  for  the  DS/CDMA  system:  (I)  perfect  channel  estimation  and  ML  data 
detection,  (II)  MMSE/ML  receiver,  and  (III)  SM/MRC  receiver.  Number  of  training  bits:  10  {dash- 
dot  lines),  30  {dashed),  50  {solid).  SIR  =  — 15  dB 


To  obtain  an  intuitive  feeling  for  what  is  taking  place,  consider  the  three  curves 
shown  in  Fig.  9.4.  All  three  correspond  to  the  use  of  the  MMSE/ML  receiver,  and 
all  three  show  the  transfer  function  of  that  receiver.  A  frequency  range  of  —250  to 
+250  MHz  is  shown,  and  there  is  a  narrowband  Gaussian  interferer  located  between 
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—20  and  +20  MHz.  The  curves  correspond  to  SIRs  of  —20  dB,  — 10  dB,  and  infinity 
(i.e.,  no  interference).  Note  that  the  three  curves  virtually  lie  on  top  of  one  another 
over  the  entire  frequency  range  shown,  except  for  the  region  —20  to  +20  MHz.  In 
that  region,  it  can  be  seen  that  the  receiver  is  implementing  a  notch  to  suppress  the 
interference  (when  it  is  present),  and  the  depth  of  that  notch  is  proportional  to  the 
strength  of  the  interference. 


Fig.  9.4  Frequency  response  of  the  ML  detector  in  the  absence  of  NBI  {solid  line)  and  in  the 
presence  of  NBI  with  SIR  =  — 10  dB  {dashed  line)  and  —20  dB  {dotted  line).  The  NBI  has  40  MHz 
of  bandwidth  and  is  centered  at  0  Hz 


9.3  Interference  Mitigation  in  Block-Modulated  Multicarrier 
Systems 

In  this  section,  two  multicarrier  schemes  [7,  15],  specifically  multicarrier  code¬ 
division  multiple  access  (MC-CDMA)  [10,  16,  23,  59]  and  orthogonal  frequency 
division  multiplexing  (OFDM)  [11,  54],  are  evaluated  in  the  presence  of  narrow- 
band  interference  (NBI).  This  interference  can  arise  when  other  systems  are  occu¬ 
pying  the  same  spectrum,  as  occurs  naturally  in  the  unlicensed  bands  [27]  or  from 
intentional  jamming.  This  has  led  to  the  emergence  of  cognitive  radio  systems  [35]. 
As  discussed  in  the  previous  section,  the  transmitter  and  receiver  of  these  systems 
must  coordinate  which  frequencies  are  unaffected  by  the  interference.  For  effective 
operation  without  cooperation,  systems  are  limited  to  dealing  with  interference  at 
the  receiver,  which  is  what  is  examined  in  this  section. 
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MC-CDMA  and  OFDM  are  designed  as  block  modulation  schemes  that  convert 
an  input  into  multiple  low-rate  streams  that  are  transmitted  on  separate  narrowband 
subcarriers,  allowing  the  subcarriers  to  overlap  orthogonally  in  frequency,  thereby 
providing  an  increase  in  spectral  efficiency.  This  design  also  allows  the  bandwidth 
of  each  subcarrier  to  approximate  the  coherence  bandwidth  of  the  channel,  ensuring 
that  each  tone  experiences  flat  fading.  This  also  simplifies  the  equalization  process 
at  the  receiver.  These  systems  fit  naturally  in  the  context  of  cognitive  radios  because 
they  have  the  ability  to  turn  off  subcarriers  that  are  affected  by  interference  [3]. 
In  MC-CDMA,  each  subcarrier  of  a  given  user  is  multiplied  by  a  single  chip  of 
a  spreading  sequence  [39].  At  the  receiver,  after  demodulation  of  the  subcarriers, 
the  chips  are  then  correlated  with  a  locally  generated  chip  sequence  to  provide  the 
multiple-access  capability.  Conversely,  OFDM  does  not  spread  the  data  symbols 
across  the  subcarriers.  OFDM  is  thus  a  special  case  of  MC-CDMA,  in  which  the 
spreading  codes  are  simply  the  columns  of  the  identity  matrix.  The  result  of  this  is 
that  uncoded  OFDM  does  not  provide  frequency  diversity. 

These  techniques  differ  from  other  wideband  systems,  such  as  direct  sequence 
CDMA  (DS-CDMA)  and  multicarrier  direct  sequence  CDMA  (MC-DS-CDMA) 
[13,28],  which  were  discussed  in  the  previous  section.  For  these  modulations  the 
bandwidth  of  the  signal  is  greater  than  the  coherence  bandwidth  of  the  chan¬ 
nel,  inducing  frequency- selective  fading.  A  rake  [39]  receiver  can  be  used  in  ad¬ 
dition  to  obtain  path  diversity,  thereby  reducing  the  deleterious  effects  of  the 
frequency- selective  multipath  fading  channel.  A  comparison  of  MC-DS-CDMA  and 
MC-CDMA  can  be  found  in  [15, 32]. 

The  inherent  frequency  diversity  of  MC-CDMA  allows  for  improved  perfor¬ 
mance  in  interference-limited,  frequency- selective  channels  [57].  Since  OFDM 
lacks  frequency  diversity,  it  must  employ  signal  processing  and  coding  techniques 
to  mitigate  the  interference.  There  has  been  a  great  deal  of  research  concerning  mit¬ 
igation  of  NBI  in  an  uncoded  OFDM  system.  Proposed  techniques  range  from  using 
orthogonal  codes  [8,21,57],  frequency-domain  cancellation  [14,37],  receiver  win¬ 
dowing  [41],  and  excision  filtering  [12].  More  recently  there  has  been  some  work 
relating  to  NBI  suppression  in  coded  OFDM  systems  [12,45,57]. 

In  this  work,  the  prediction-error  filter  (PEF)  [30, 60]  is  selected  for  interference 
mitigation  in  OFDM.  This  is  an  example  of  excision  filtering  first  described  in  [12]. 
First,  this  system  is  compared  with  MC-CDMA  in  an  uncoded  scenario.  Then  the 
addition  of  forward  error  correction  coding  for  each  system  will  be  examined.  In 
this  case,  OFDM  will  now  also  possess  frequency  diversity  due  to  the  coding  across 
the  subcarriers. 


9.3.1  Interference  Mitigation  in  an  Uncoded  Multicarrier  System 


A  block  diagram  of  the  system  model  considered  in  this  section  is  shown  in  Fig.  9.5. 
The  code  (spreading)  matrix,  B,  is  an  A  x  A  unitary  matrix  whose  columns  contain 
the  A  spreading  sequences.  The  kth  column  is  responsible  for  spreading  the  kth 
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data  symbol,  Note  that  when  B  is  the  identity  matrix  (i.e.,  B  =  I),  the  system 
becomes  identical  to  OFDM. 


n  X 


Fig.  9.5  Block  diagram  of  uncoded  MC-CDMA 


The  received  signal  is  given  by 

Y' =  F  I^HF^BD  +  X  +  j;  j  ,  (9.1) 

where  F  is  the  N  xN  DFT  matrix,  B  is  the  N  xN  code  matrix,  D  is  the  x  1  vector 
of  data  symbols  that  is  drawn  from  an  arbitrary  QAM  constellation,  x  is  an  A  x  1 
vector  of  interference  samples,  /y  is  an  A  x  1  vector  of  Gaussian  noise  samples,  and 
(•)^  is  the  Hermitian  (conjugate  transpose)  operator.  Finally,  H  is  the  A  x  A  channel 
matrix  whose  components  are  given  by  the  taps  of  a  multipath  (frequency-selective) 
channel,  as  defined  by  the  impulse  response 


L-l 

h(t)  =  J2hiHt-lTs),  (9.2) 

1=0 

where  L  is  the  number  of  multipath  components  and  Ts  is  the  original  symbol 
period.  The  channel  coefficients,  hi,  are  i.i.d.  zero-mean,  circularly  Gaussian  ran¬ 
dom  variables  with  variance 


E[\hlf]  =  l 

where  E  [•]  represents  the  expectation  operator. 

A  cyclic  prefix,  composed  of  the  last  G  samples  of  F^D,  is  employed  to  mit¬ 
igate  the  effects  of  intersymbol  interference  and  intercarrier  interference  (ICI). 
When  the  length  of  the  cyclic  prefix  is  greater  than  the  delay  spread  of  the  chan¬ 
nel  (L),  linear  convolution  is  equivalent  to  circular  convolution  and  H  is  circu- 
lant.  A  well-known  property  of  circulant  matrices  is  that  their  eigenvectors  are 
simply  the  columns  of  the  DFT  matrix  [22].  Thus,  (9.1)  can  be  rewritten  as 
follows: 


Y'  =  HBD  +  Fx  +  Fjy 

=  HBD  +  X  +  N,  (9.4) 

where  H  is  an  A  x  A  diagonal  matrix  of  the  eigenvalues  of  H.  Coherent  demodu¬ 
lation  requires  removing  the  phase  of  the  channel,  thus  let  be  a  diagonal  matrix 
containing  the  phase  information  of  the  channel,  that  is 
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Hfl 


jZh 


(9.5) 


and  let  Hp  be  the  magnitude  of  the  eigenvalues,  given  as 


|H| 


po  0  O' 

0  Pi  •••  0 

0  0  •  • • pn-i 


where  pk  is  a  Rayleigh  random  variable.  Note  that  H  =  H@Hp. 

After  coherent  detection,  the  input  to  the  combiner  is  given  as 

Y  =  Hf  Y'  =  HpBD  +  Hf  (X  +  N) . 

The  interference  samples  in  the  time  domain  seen  in  (9.1)  are  given  by 

Xn  = 


(9.6) 


(9.7) 


(9.8) 


where  Ex  is  the  interferer  power  and  ^  is  a  random  phase  distributed  uniformly 
between  — tt  and  tt.  The  frequency  of  the  interference,  fx,  is  defined  as  a  function 
of  the  original  sampling  frequency,  fs,  namely  fx  =  (m  +  a)  In  this  formulation 
of  the  interference,  m  represents  the  tone  closest  to  the  interference,  while  a  G 
[—0.5,  0.5]  is  the  offset  from  tone  m.  The  second  line  of  (9.8)  is  obtained  by  noting 
that  fs  =  Y,- 

The  time-domain  interference  samples  in  (9.8)  are  converted  into  the  frequency 
domain 


Xk  =  [Fx]^ 


V  N  I  _ 


(9.9) 


where  is  the  ^th  component  of  the  given  vector.  The  result  of  (9.9)  indicates 
that  the  interference  per  tone  is  dependent  on  the  relative  distance  of  tone  k  from 
the  interferer  frequency.  This  is  referred  to  as  spectral  leakage,  as  the  interference 
power  is  leaked  into  all  of  the  tones.  From  (9.9),  it  can  be  seen  that  as  the  interference 
power  increases,  so  does  the  number  of  tones  that  are  significantly  interfered  with. 

In  the  case  of  a  =  0,  the  interference  is  orthogonal  to  all  the  subcarriers  other 
than  the  subcarrier.  The  interference  after  the  DFT  reduces  to 


Xk 


k  =  m 

0,  k  ^  m 


(9.10) 


Thus  the  interference  is  limited  to  tone  m  in  this  case. 

Finally,  the  last  term  of  (9.4),  N,  is  an  x  1  vector  of  noise  samples  that  have 
the  same  statistics  as  in  the  time  domain. 
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9.3.1.1  MC-CDMA 

The  code  matrix  (B)  employed  in  this  work  is  defined  component  wise  [36, 57]  by 

[B]h  =  k,l  =  0,...,N  -1.  (9.11) 

Note  that  the  code  matrix  is  unitary,  that  is  BB^  =  B^B  =  In,  where  In  is  the 
N  xN  identity  matrix.  Another  option  that  can  be  used  for  the  orthogonal  spreading 
codes,  are  the  Walsh-Hadamard  codes.  However,  note  that  these  codes  are  only 
defined  for  lengths  of  A  =  2^,  where  n  is  an  integer  greater  than  zero  [15]. 

The  combiner  weights  required  to  extract  the  symbol  vector,  D,  are  derived  under 
the  minimum  mean-square  error  (MMSE)  criterion  [57,58]  and  obtained  by  solving 
the  Wiener-Hopf  equation 

W  =  (9.12) 

where  Ryy  is  the  autocorrelation  matrix  of  the  received  signal  and  Ryd  is  the  cross¬ 
correlation  matrix  between  the  received  signal  and  the  symbol  vector,  D.  Assuming 
that  all  the  received  components  are  independent  and  that  the  channel  amplitudes 
(pk)  are  known  and  deterministic,  the  autocorrelation  matrix  is  given  by 


Ryy  =  Ebiipiip  +  NqIn  +  Rx,  (9.13) 


where  Rx  is  the  autocorrelation  matrix  of  the  interference  samples  in  the  frequency 
domain  defined  as 

Rx  =HfFR;,F"H0,  (9.14) 

where  Rx  is  the  time-domain  correlation  matrix  of  the  interference,  its  entries 
given  by 

rAD  =  (9.15) 

Similarly,  the  cross-correlation  vector  is  given  by 

Ryd  =  EtiipB.  (9.16) 


For  the  case  of  no  interference  (Rx  =  0),  the  weights  are  given  as 


W  =  diag 


PoEb 


PlEb 


P^Eb  +  Aq  P? Eb  +  Aq 


PN-lEb 


|B, 


(9.17) 


where  diag(-)  returns  a  diagonal  matrix  from  the  inputted  vector. 

When  the  interference  is  orthogonal  to  the  subcarriers  (except  for  the  sub¬ 
carrier),  O'  =  0,  recall  that  X  will  have  one  non-zero  component,  at  index  m.  Thus, 
Rx  will  also  only  have  one  non-zero  value,  equal  to  N Ex  at  position  (m,  m).  Let 
w  =  [w(0),  . . . ,  w(N  —  1)],  where  w(l)  is  given  by 
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w(l)  = 


PlEb 

pfEb+No  ’ 

PlEb 

pf  Eb+No+N  Ex  ’ 


I  ^  m 
I  =  m 


The  MMSE  weights  are  then  given  as 


(9.18) 


W  =  diag(w)B.  (9.19) 

It  can  be  seen  from  (9.19)  that  when  interference  is  not  present  on  the  tone,  the 
weight  is  the  same  as  for  the  no-interference  case.  Note  that  when  I  =  m  and  the 
narrowband  interference  is  strong,  then  Wk(l)  is  essentially  zero,  thus  ignoring  the 
subcarrier  that  is  unreliable. 

Finally,  for  the  case  of  a  nonorthogonal  interfere^  the  components  of  the  auto¬ 
correlation  matrix  for  the  interference  are  given  by 


Ex  2  —  2cos(27rQ')  ^  ^  ^ 

N  I  _  ^j^im+a-k)  _  ^-j^{m+a-l)  ^j^Q-k)  ’ 


(9.20) 

To  find  the  weights  one  could  use  (9.20)  and  invert  (9.12).  However,  noticing  that 
Rx  is  rank-1,  having  one  dominant  eigenvalue,  then  Rx  can  be  approximated  using 
the  singular  value  decomposition  (SVD)  [25]  as  follows: 


Rx  =  USV" 


(9.21) 


where  X  =  N Ex  is  the  non-zero  eigenvalue  of  Rx,  u  is  the  first  column  of  U,  and  v 
is  the  first  column  of  V.  Then  using  the  Sherman-Morrison  formula  [25],  the  MMSE 
combiner  weights  are  derived  to  be 


W  = 


AAuv^A 
1  +  kv^Au 


EbUpB, 


(9.22) 


where  A  =  diag 


p^Eb-\-No  pfEb+No  pN_iEb+No 

Finally,  the  estimates  of  the  transmitted  data  symbols  are  given  by 


D 


W^Y. 


(9.23) 


Note  that  a  consequence  of  MMSE  combining  in  frequency-selective  channels  is 
that  there  is  a  loss  of  orthogonality  among  the  data  symbols, 

D  =  D  +  -  lx)  D  +  (X  +  N) .  (9.24) 

On  the  other  hand,  when  the  channel  is  AWGN  or  fiat  fading,  equal  gain  combining 
(i.e.,  W  =  B)  can  be  utilized  to  preserve  orthogonality,  that  is 

D  =  D  +  (X  +  N) . 


(9.25) 
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9.3.1.2  OFDM 

Unlike  MC-CDMA,  OFDM  does  not  provide  frequency  diversity  since  the  code 
matrix  in  Fig.  9.5  is  set  to  the  identity  matrix.  The  prediction-error  filter  (PEF)  is 
considered  for  this  system  as  a  means  for  removing  the  interference  in  the  time 
domain,  thereby  avoiding  the  spectral  leakage  that  occurs  after  demodulation  (see 
block  diagram  given  in  Fig.  9.6).  The  PEF  is  a  well- studied  structure  that  uses  the 
correlation  between  past  samples  to  form  an  estimate  of  the  current  sample  [40]. 
The  PEF  has  the  property  that  it  removes  the  correlation  between  samples,  thereby 
whitening  the  spectrum. 


V  X 


Fig.  9.6  Block  diagram  of  uncoded  OFDM  with  the  PEF 

This  technique  has  been  used  to  remove  narrowband  interference  in  many  appli¬ 
cations,  such  as  spread-spectrum  systems  when  the  processing  gain  does  not  provide 
enough  immunity  to  the  interference  [49, 53].  A  block  diagram  of  the  PEF  is  shown 
in  Fig.  9.7.  It  is  simply  a  transversal  filter  with  M  taps.  The  decorrelation  delay  (A) 
ensures  that  the  current  sample  is  decorrelated  from  the  samples  in  the  filter  when 
calculating  the  error  term. 


Fig.  9.7  Block  diagram  of  the  prediction-error  filter  (PEF) 


The  output  of  the  PEF,  is  defined  as  the  subtraction  of  the  estimate  of  the 
interference  from  the  current  input  sample. 


M-l 

Zn=yn-  cfyn-A-h  (9.26) 

1=0 

where  c/  are  the  weights  of  the  predictor  and  (•)*  represents  complex  conjugation. 
The  optimal  weights  under  the  minimum  mean- squared  error  (MMSE)  criterion  are 


230 


Amn  Batra,  James  R.  Zeidler,  John  G.  Proakis,  and  Laurence  B.  Milstein 


found  using  the  orthogonality  principle  and  the  method  of  undetermined  coefficients 
[61].  The  weights  are  given  by 

C,  = /  =  0,  (9.27) 

where  K  is  given  by 

Ex 

K  = - ^ - .  (9.28) 

+  yVo  +  me. 

For  the  scenario  of  interest  in  this  chapter,  the  interference  power  is  much  larger 
than  both  the  signal  power  and  the  noise  power.  Therefore,  the  SIR  and  the  noise  to 
interference  (NIR)  can  be  assumed  to  be  very  small  (i.e.,  SIR  ^  0  dB,  NIR  ^  0  dB 
[31])  allowing  the  coefficient,  to  be  approximated  as  ^  ~ 

This  filter  is  easily  implemented  adaptively  using  the  least  mean  square  (LMS) 
algorithm.  The  convergence  properties  are  well  described  [5,42].  Note  also  that  Zn 
is  the  error  term  for  the  adaptive  structure,  indicating  that  no  training  symbols  are 
required  when  calculating  the  error  term. 

The  PEF  is  implemented  before  the  removal  of  the  cyclic  prefix.  Let  the  one- step 
predictor  weights  be  defined  as  c  =  [l,  —  Cq,  •  •  • .  and  the  convolution  of 

the  filter  and  channel  be  defined  as  a  =  c  *  h.  It  is  assumed  that  the  overall  length 
of  a  is  less  than  the  length  of  the  cyclic  prefix  (i.e.,  L  -h  M  —  1  <  G)  to  ensure  that 
there  is  no  ISI  or  ICI.  Therefore,  the  effective  channel  matrix.  A,  is  circulant.  Then, 
the  filtered  signal  in  the  frequency  domain  can  be  written  as 

Z  =  FAF  D  +  FCnoise  (^cp  T  ^cp) 

=  ad  +  FCnoise  (^cp  H”  ^cp)  ’  (9.29) 

where  A  is  the  diagonal  matrix  of  the  eigenvalues  of  A,  Xcp  and  are  vectors  of 
interference  and  noise  samples  (length  A  +  G),  respectively  that  are  not  cyclically 
extended,  and  Cnoise  is  the  A  x  (A+G)  filtering  matrix  for  the  noise  and  interference 
that  is  defined  as 

Cnoise  =  G-M  Toeplitz  [c*,  j  ,  (9.30) 

where  0/,^  is  the  i  x  j  zero  matrix.  The  Toeplitz  operator,  Toeplitz (c6>/Mm^,  row), 
generates  a  Toeplitz  matrix  from  a  column  vector  and  a  row  vector.  Note  that  A  in 
(9.32)  can  also  be  defined  as 

A  =  HC,  (9.31) 

where  C  =  ^/N¥  [c*  0i,a^_(m+1)]^.  This  vector  is  the  sampled  frequency  response 
of  the  notch  filter,  with  the  notch  located  on  the  tone  closest  to  the  interference. 

For  coherent  demodulation,  the  phase  of  the  effective  channel  matrix  A  must  be 
removed;  therefore,  similar  to  (9.5)  and  (9.6),  let  A  =  AgAp,  where  Ap  =  |A|  and 

Aq  =  Thus  for  coherent  demodulation. 
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z  =  Af  Z'  =  ApD  +  Af  FCnoise  (xcp  +  »/cp)  •  (9.32) 

Note  that  the  diagonal  elements  of  Ap  are  not  necessarily  distributed  as  Rayleigh 
random  variables  due  to  the  PER 

Let  the  uncanceled  interference  and  noise  be  grouped  into  one  general  noise  term 

N  =  Ap  FCnoise  (^cp  +  ^cp)  *  (9.33) 

It  is  clear  from  (9.33)  that  the  noise  samples,  Nk^  are  correlated  due  to  the  PEE 
matrix,  Cnoise-  H  is  also  noted  that  the  noise  samples  Nk  are  assumed  to  be  Gaussian 
random  variables.  As  stated  in  [31],  this  system  is  difficult  to  analyze  when  the 
noise  samples  are  not  strictly  independent,  however,  due  to  the  fact  that  the  noise 
power  in  the  main  tap  (equal  to  unity)  is  much  larger  than  in  the  remaining  taps 
(approximately  ^  for  the  scenario  of  interest  in  this  chapter),  it  is  reasonable  to 

assume  that  the  noise  samples  are  independent.  Therefore,  let  be  the  variance 
of  the  noise  on  tone  k,  given  by 


Note  that  the  variances  given  in  (9.34)  are  scaled  according  to  the  notch  filter  that  is 
used  to  suppress  the  interference. 

One-tap  equalization  is  performed  instead  of  MMSE  combining  as  is  done  for 
MC-CDMA.  Thus  the  estimates  of  the  transmitted  data  symbols  are  given  by 

D  =  Z/Ap.  (9.35) 

An  analytic  expression  for  the  BER  of  this  OEDM  system  given  knowledge  of  the 
channel  amplitudes  is  simply  the  average  of  the  per-tone  BER,  given  by 


Pe 


1 

N 


N-l 

Eq 


k=0 


CT?  (k) 


(9.36) 


The  analytical  results  of  (9.36)  are  seen  in  Eig.  9.8  along  with  the  simulation  results 
for  SIR  =  — 20dB  and  N  =  64,  G  =  16,  L  =  5,  and  M  =  12.  The  simulation  and 
theoretical  results  agree  quite  well,  validating  the  assumption  of  independent  noise 
samples.  Also  note  that  an  error  floor  arises,  because  the  OEDM  system  is  limited  by 
its  worst  set  of  tones.  In  this  case,  these  tones  are  the  ones  located  around  the  notch 
generated  by  the  PEE.  This  suggests  the  need  for  forward  error  correction  coding. 
The  addition  of  a  convolutional  code  can  reduce  the  BER  from  10“^  to  10“^.  Note 
that  without  the  use  of  the  PEE,  the  BER  is  so  poor  that  forward  error  correction 
coding  cannot  be  effectively  employed. 
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Fig.  9.8  Uncoded  BER  results  for  theoretical  and  simulated  OFDM  with  the  prediction-error  filter 
(PEE).  Also  plotted  are  the  theoretical  BER  curves  for  AWGN  and  Rayleigh  fading.  Parameters: 
=  64,  G  =  16,  L  =  5,  M  =  12,  SIR  =  -20dB 


9.3.1.3  Results 

A  comparison  of  the  two  techniques  described  previously  is  compared  in  an  uncoded 
scenario  in  Fig.  9.9.  For  reference,  the  theoretical  BER  curves  for  both  AWGN  and 
Rayleigh  fading  are  provided.  Each  system  is  simulated  for  the  cases  (i)  no  inter¬ 
ference  and  (ii)  one  nonorthogonal  interferes  The  plot  for  MC-CDMA  in  no  inter¬ 
ference  indicates  that  frequency  diversity  provides  a  benefit  in  a  frequency- selective 
channel.  The  performance  of  MC-CDMA  with  a  single  nonorthogonal  interferer  is 
very  close  to  the  case  of  no  interference.  In  this  case,  the  combiner  utilizes  the  statis¬ 
tics  of  the  channel  and  interference  to  successfully  recombine  the  transmitted  signal. 
The  deviation  from  the  case  of  no  interference  is  very  small. 

The  PEE  used  with  OFDM  assumes  knowledge  of  the  MMSE  weights  given  in 
(9.27).  In  the  presence  of  no  interference  the  performance  of  the  system  reduces  to 
that  of  Rayleigh  fading,  as  is  expected.  When  interference  is  present,  the  case  of  a 
nonorthogonal  interference  approximates  the  case  of  Rayleigh  fading  at  low  Ejj  /  Aq 
values;  however  at  high  Eij/Nq  values,  an  error  floor  arises.  This  is  due  to  the  notch 
filter  that  removes  a  few  tones  when  mitigating  the  interference. 

For  both  MC-CDMA  and  OFDM  with  PEF,  the  interference  must  be  estimated 
in  order  to  be  implemented  in  a  real-world  communication  system.  As  mentioned 
earlier,  the  PEF  is  easily  implemented  using  the  low-complexity  LMS  algorithm 
without  training  symbols.  This  allows  the  interference  to  be  adaptively  estimated 
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Fig.  9.9  Uncoded  BER  results  for  MC-CDMA  and  OFDM  with  the  prediction-error  filter  (PEE) 
for  the  cases  of  (i)  no  interference  and  (ii)  one  nonorthogonal  interferes  Also  plotted  are  the 
theoretical  BER  curves  for  AWGN  and  Rayleigh  fading.  Parameters:  N  =  64,  G  =  16,  L  =  5, 
M  =  12,  SIR=  -20  dB 


and  removed.  In  the  case  of  MC-CDMA,  the  correlation  matrix  of  the  interference  is 
required  when  determining  the  MMSE  combiner  weights.  This  can  be  accomplished 
using  training  symbols  to  provide  an  estimate  for  the  received  signal  correlation 
matrix. 


9.3.2  Interference  Mitigation  in  Coded  Multicarrier  Systems 


In  this  section,  the  impact  of  applying  forward  error  correction  to  the  uncoded 
systems  analyzed  above  is  evaluated.  Block  diagrams  for  coded  MC-CDMA  and 
coded  OFDM  with  PEF  are  depicted  in  Figs.  9.10  and  9.11,  respectively.  A  con¬ 
volutional  code  of  rate  1  /2,  constraint  length  equal  to  7  with  generating  polynomi¬ 
als  (133,  171)8  is  used.  One  codeword  spans  50  MC-CDMA/OFDM  symbols.  At 
the  receiver  the  Viterbi  [39]  decoder  is  used  to  decode  the  codeword  utilizing  soft 
decisions.  Note  that  a  block  interleaver  is  used  with  OFDM  to  provide  frequency 
diversity.  An  interleaver  can  also  be  used  with  MC-CDMA  when  the  data  are  spread 
over  a  subset  of  the  subcarriers. 

The  soft  decisions  inputted  into  the  Viterbi  decoder  encourage  the  use  of  erasures 
when  reliable  log-likelihood  ratios  (FFRs)  are  unavailable.  In  OFDM,  erasures  are 
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Fig.  9.10  Block  diagram  for  coded  MC-CDMA 


Fig.  9.11  Block  diagram  for  coded  OFDM  with  the  PEF 


only  effective  when  the  interference  is  weak  [4].  However,  when  the  interference  is 
strong,  the  spectral  leakage  of  interference  power  causes  a  large  number  of  tones  to 
be  corrupted.  Insertion  of  erasures  ends  up  compromising  the  code’s  error  correc¬ 
tion  capability.  The  use  of  the  PEF  can  be  considered  as  a  form  of  erasure  insertion, 
whereby  the  erasures  are  localized  around  the  interference  by  the  notch  that  is  placed 
in  the  frequency  spectrum  of  the  received  signal.  The  code  is  then  tasked  with  cor¬ 
recting  any  erasures  caused  by  the  channel  or  interference.  The  full  erasures  are 
inserted  after  the  PEF  using  the  Bayesian  erasure  formula  described  in  [6]. 

The  LLR  values  [15]  that  are  inputted  into  the  Viterbi  decoder  for  MC-CDMA 
with  BPSK  modulation  are  given  as 


^MC-CDMAiDk\Yk)  = 


2mYk} 


(9.37) 


where  94{-}  is  the  real  operator  and  is  the  variance  of  the  ICI,  noise,  and  inter¬ 
ference  as  defined  in  (9.24).  Assuming  these  vectors  are  composed  of  independent 
Gaussian  random  variables,  is  given  as  follows: 


(7^  =  diag 


Eh  (w^HpB  - 


+  W"H^FR;,F"HeW 


(9.38) 


Similarly,  for  OFDM 


2ioFDM(^y^l^y^)  = 


2mZk} 


<7^ 


(9.39) 


where  is  given  as 


diag  FCnoise 

.  +  f  >»)  C,*!"'*.] 

diag 

1  1 

1 _ 1 

(9.40) 
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9.3.2. 1  Results 

Figure  9.12  demonstrates  the  results  for  the  coded  simulation  of  MC-CDMA  and 
OFDM  with  the  PEF  for  the  cases  of  (i)  no  interference  and  (ii)  one  nonorthogonal 
interferer,  for  A  =  64,  G  =  16,  L  =  5,  M  =  12,  SIR  =  — 20dB.  Also  plotted 
are  the  theoretical  BER  curves  for  uncoded  Rayleigh  fading  and  a  BER  bound  for 
AWGN.  It  can  be  clearly  seen  that  the  performance  of  both  systems  with  interference 
is  very  close  to  the  case  of  no  interference.  Also  note  that  the  performance  of  OFDM 
with  frequency  diversity  provided  by  coding  and  the  PEF  as  an  erasure  insertion 
technique  is  equivalent  to  that  of  coded  MC-CDMA  for  which  frequency  diversity 
is  obtained  through  spreading.  The  coding  gain  is  also  apparent  when  comparing 
the  coded  results  with  Fig.  9.9. 


Fig.  9.12  Coded  BER  results  for  MC-CDMA  and  OFDM  with  the  prediction-error  filter  (PEF)  for 
the  cases  of  (i)  no  interference  and  (ii)  one  nonorthogonal  interferer.  Also  plotted  are  the  theoretical 
BER  curves,  uncoded  Rayleigh  Fading,  and  a  BER  bound  for  AWGN.  Parameters:  N  =  6A,  G  = 
16,  L  =  5,  M  =  12,  SIR  =  -20  dB 


9.33  Doppler  Sensitivity  of  OFDM  in  Mobile  Applications 

Mobile  systems  that  employ  OFDM  may  suffer  degradations  from  ISI  and  ICI  in 
addition  to  the  other  forms  of  interference  discussed  above.  In  [50],  the  Doppler  sen¬ 
sitivity  of  OFDM  was  investigated  for  PSK  and  QAM  modulations  in  the  presence 
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of  ICI  and  the  loss  in  performance  due  to  Doppler  spreading  was  evaluated.  Also,  in 
[52],  comparisons  were  made  between  the  Doppler  sensitivity  of  OFDM  and  filtered 
multitone  (FMT)  modulation  in  frequency-selective  time- varying  channels. 

FMT  is  an  alternate  form  of  multicarrier  modulation  that  has  been  proposed 
for  very  high-rate  subscriber  lines  and  for  wireless  communications.  FMT  avoids 
spectral  overlap  between  subcarriers  by  implementing  a  noncritical  sampling  tech¬ 
nique  that  satisfies  a  perfect  reconstruction  condition,  so  that  the  sampling  signals 
at  the  receiver  are  free  from  ISI  as  well  as  ICI.  The  ISI  and  ICI  in  an  FMT  sys¬ 
tem  is  suppressed  at  the  expense  of  choosing  a  larger  upsampling-factor-to-number- 
of- subcarriers  ratio.  Due  to  the  phase  and  amplitude  distortion  introduced  by  the 
fading  channel,  not  only  is  the  orthogonality  among  different  subcarriers  destroyed 
but  also  the  perfect  Nyquist  sampling  condition  of  the  baseband  matched  filters 
is  no  longer  valid.  Consequently,  ICI  as  well  as  ISI  will  cause  distortions  to  the 
transmitted  signal.  The  interference  caused  by  the  channel  frequency  selectivity 
and  time  variance  was  quantified  by  analyzing  the  demodulated  signals  at  the  re¬ 
ceiver  under  different  fading  channel  conditions.  It  was  shown  [52]  analytically  that 
quasi- static  frequency  fiat  fading  channels  do  not  introduce  ICI  or  ISI  to  a  multi¬ 
carrier  system,  while  frequency- selective  channels  introduce  both  ISI  and  ICI  in  an 
FMT  system.  The  resulting  carrier-to-interference  (C/I)  ratio  obtained  under  differ¬ 
ent  channel  conditions  was  shown  to  provide  a  trade-off  between  spectral  efficiency 
and  system  performance  degradation.  Also  it  was  shown  [52]  that  OFDM  outper¬ 
forms  FMT  in  a  static  or  slowly  fading  and  low  frequency- selective  channel  but  is 
inferior  to  FMT  modulation  in  highly  frequency- selective  and  time- varying  fading 
environments. 

This  work  was  extended  in  [51]  to  consider  a  per-channel-equalized  FMT  ar¬ 
chitecture.  It  is  shown  that  with  a  sufficient  number  of  equalizer  coefficients  it  is 
possible  to  mitigate  the  ISI,  but  the  ICI  remains.  A  closed- form  analytical  expres¬ 
sion  for  the  C/I  ratio  and  its  upper  bound  is  derived  and  the  trade-offs  in  spectral 
efficiency  and  system  performance  are  further  quantified. 


9.4  Interference  Suppression  in  Broadcast  MIMO  Systems 

In  this  section,  we  consider  multiple-input,  multiple-output  (MIMO)  wireless  com¬ 
munication  systems  that  employ  multiple-transmit  and  multiple-receive  antennas  to 
increase  the  data  rate  and  achieve  signal  diversity  in  fading  multipath  channels.  The 
performance  of  MIMO  systems  is  degraded  by  two  types  of  interference.  One  is 
intersymbol  interference  due  to  channel  multipath  dispersion.  The  other  is  cross¬ 
talk  or  interchannel  interference  due  to  the  simultaneous  transmissions  from  the 
multiple-transmit  antennas.  The  focus  of  this  section  is  on  equalization  and  detection 
algorithms  for  mitigating  these  two  types  of  interference.  In  particular,  we  consider 
a  point-to-multipoint  (broadcast)  MIMO  transmission  system  in  which  the  chan¬ 
nel  characteristics  are  assumed  to  be  known  at  the  transmitter,  so  that  interference 
mitigation  is  performed  at  the  transmitter. 
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Let  us  consider  a  point-to-multipoint  MIMO  system  which  transmits  data  simul¬ 
taneously  to  multiple  users  that  are  geographically  distributed.  The  transmitter  is 
assumed  to  employ  Nt  antennas  to  transmit  data  to  K  receivers,  where  Nt  >  K. 
Each  user  is  assumed  to  have  a  receiver  with  one  or  more  receiving  antennas.  This 
scenario  applies,  for  example,  to  the  downlink  of  a  wireless  local  area  network.  The 
distinguishing  feature  of  the  MIMO  broadcast  system  is  that  the  receivers,  which 
are  geographically  separated,  do  not  employ  any  coordination  in  processing  the  re¬ 
ceived  signals. 

In  a  MIMO  broadcast  system,  there  are  two  possible  approaches  for  dealing  with 
the  multiple-access  interference  (MAI)  resulting  from  the  simultaneous  transmis¬ 
sion  to  multiple  users.  One  approach  is  to  have  each  receiver  employ  interference 
mitigation  in  the  recovery  of  its  desired  signal.  In  most  cases,  this  approach  is  im¬ 
practical  because  the  users  lack  the  processing  capability  and  are  constrained  by 
the  limited  energy  resources  inherent  in  the  use  of  battery  power.  The  alternative 
approach  is  to  employ  interference  mitigation  at  the  transmitter,  which  possess  sig¬ 
nificantly  more  processing  capability  and  energy  resources. 

To  mitigate  the  MAI  at  the  transmitter,  the  transmitter  must  know  the  channel 
characteristics,  typically  the  channel  impulse  response.  This  channel  state  infor¬ 
mation  (CSI)  may  be  obtained  from  channel  measurements  performed  at  each  of 
the  receivers  by  means  of  received  pilot  signals  sent  by  the  transmitter.  Then,  the 
CSI  must  be  sent  to  the  transmitter.  In  such  a  scenario,  the  channel  time  variations 
must  be  relatively  slow  so  that  a  reliable  estimate  of  the  channel  characteristics  is 
available  at  the  transmitter.  In  some  systems,  the  uplink  and  downlink  channels  are 
identical,  e.g.,  the  same  frequency  band  is  employed  for  both  the  uplink  and  the 
downlink,  but  separate  time  slots  are  used  for  transmission.  This  transmission  mode 
is  called  time-division  duplex  (TDD).  In  TDD  systems,  the  pilot  signals  for  channel 
measurement  may  be  sent  by  each  of  the  users  in  the  uplink.  In  our  treatment  below, 
we  assume  that  the  CSI  at  the  transmitter  is  perfect. 

The  suppression  of  MAI  by  means  of  transmitter  processing  is  usually  called 
signal  precoding  or  signal  preprocessing.  Signal  precoding  at  the  transmitter  may 
take  one  of  several  forms,  depending  on  the  criterion  or  the  method  used  to  perform 
the  precoding.  The  simplest  precoding  methods  are  linear  and  are  based  on  either 
the  zero-forcing  (ZF)  criterion  or  the  MSE  criterion.  On  the  other  hand,  there  are 
nonlinear  signal  precoding  methods  that  result  in  better  system  performance.  First, 
we  treat  linear  precoding  and,  then,  we  describe  nonlinear  precoding  methods  and 
compare  their  performance. 


9.4.1  Linear  Precoding  of  the  Transmitted  Signals 

In  this  section,  for  mathematical  convenience,  we  assume  that  each  user  has  a  single 
antenna  and  the  number  of  receivers  (users)  is  ^  <  Nj.  It  is  also  convenient  to 
assume  that  the  channel  is  nondispersive.  The  communication  system  configuration 
is  shown  in  Fig.  9.13,  where  the  precoding  matrix  is  denoted  as  A^.  Hence,  the 
received  signal  vector  is 
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y  =  HA7S  +  jy,  (9.41) 

where  H  is  an  ^  x  Nt  matrix,  At  is  an  Nt  x  K  precoding  matrix,  s  is  a  ^  x  1  data 
vector,  and  jy  is  a  ^  x  1  Gaussian  noise  vector.  The  matrix  that  eliminates  the  MAI 
at  each  receiver  is  generally  given  by  the  Moore-Penrose  pseudo-inverse 

H+ =H"(HH^)“*.  (9.42) 


Hence,  the  precoding  matrix  is 


At-  =  kH+,  (9.43) 

where  y  is  sl  scale  factor  that  is  selected  to  satisfy  the  total  transmitted  power 
allocation,  i.e.,  ||  A^s  |p=  P.  Thus,  the  precoding  matrix  in  (9.43)  allows  the  in¬ 
dividual  users  to  recover  their  desired  symbols  without  any  interference  from  the 
signals  transmitted  to  the  other  users.  We  also  observe  that  in  the  special  case  where 
K  =  Nt,  A^  =  so  that  the  precoding  matrix  is  proportional  to  the  in¬ 

verse  channel  matrix.  This  constitutes  a  zero-forcing  equalizer  implemented  at  the 
transmitter. 


Fig.  9.13  Model  of  MIMO  broadcast  system  employing  linear  precoding 


Figure  9.14  illustrates  the  error  rate  performance  of  the  zero-forcing  precoder 
obtained  via  Monte  Carlo  simulation  for  K  =  Nt  =  ^,  6,  and  10  and  QPSK 
modulation.  The  channel  matrix  elements  are  realizations  of  complex-valued  i.i.d. 
zero-mean  Gaussian  random  variables  with  unit  variance.  We  observe  that  the  error 
rate  increases  with  an  increase  in  the  number  of  users.  We  attribute  this  deteriora¬ 
tion  in  performance  to  the  ill-conditioning  of  the  channel  matrix  H.  This  is  the  major 
drawback  with  the  zero-forcing  precoder. 

If  we  relax  the  condition  that  the  interference  be  zero  at  all  the  receivers,  the 
performance  degradation  can  be  reduced.  This  can  be  accomplished  by  using  the 
linear  MSB  criterion  in  the  design  of  the  precoding  matrix  A^.  Thus,  we  select  A^ 
to  minimize  the  cost  function 


J(At,  y)  =  arg  min  E 
yNt 


—  (HA^^s  +  jy)  —  s 

y 


(9.44) 
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Fig.  9.14  Performance  of  ZF  linear  precoding  with  Nj  =  K  =  4,  6,  10 


subject  to  the  transmitted  power  allocation  ||Ars||^  =  P  and  where  the  expecta¬ 
tion  in  (9.44)  is  taken  over  the  noise  and  signal  statistics.  The  solution  to  the  MSE 
criterion  is  the  precoding  matrix 

Ar  =  ]/H^(HH^  -h  pi)-\  (9.45) 

where  y  is  the  scale  factor  that  is  selected  to  satisfy  the  power  allocation  and  ^ 
is  defined  as  a  loading  factor,  which  when  selected  as  ^0  =  K/P  maximizes  the 
signal-to-interference-plus-noise  ratio  (SINR)  at  the  receiver  [38]. 

The  error  rate  performance  of  the  MMSE  linear  precoder  obtained  by  Monte 
Carlo  simulation  in  a  frequency-nonselective  Rayleigh  fading  channel  is  illustrated 
in  Eig.  9.15  for  K  =  Nt  =  ^,6,  and  10.  We  observe  that  the  error  rate  perfor¬ 
mance  improves  slightly  as  the  number  of  users  K  increases  and  that  it  exceeds  the 
performance  of  the  zero-forcing  precoder. 


9.4.2  Nonlinear  Precoding  of  the  Transmitted  Signals:  The  QR 
Decomposition 


When  the  transmitter  knows  the  interference  caused  on  other  users  by  the  trans¬ 
mission  of  a  signal  to  any  particular  user,  the  transmitter  can  design  signals  for 
each  of  the  other  users  to  cancel  the  interference.  The  major  problem  with  such  an 
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Fig.  9.15  Performance  of  MMSE  linear  precoding  with  Nt  =  K  =  4,  6,  10 


approach  is  to  perform  the  interference  cancellation  without  increasing  the  transmit¬ 
ter  power.  This  same  problem  is  encountered  in  decision-feedback  channel  equal¬ 
ization  of  a  single  user,  where  the  feedback  filter  of  the  DFE  is  implemented  at 
the  transmitter.  In  that  case,  when  the  difference  between  the  desired  transmit¬ 
ted  symbol  and  the  ISI  exceeds  the  range  of  the  desired  transmitted  symbol,  the 
difference  is  reduced  by  subtracting  an  integer  multiple  of  2M  for  M-ary  PAM, 
where  [— M,  M)  is  the  range  of  the  desired  transmitted  symbol.  This  same  non¬ 
linear  precoding  method  called  Tomlinson-Harashima  precoding  [24,  47]  can  be 
applied  to  the  cancellation  of  the  MAI  in  a  MIMO  broadcast  communication  system 
[17,18,56]. 

Figure  9.16  illustrates  the  precoding  operations  for  the  MIMO  multiuser  system. 
The  channel  impulse  response  between  the  ith  transmit  antenna  and  the  receive 
antenna  of  the  ^th  user  is  given  by 


L-l 

hki{t)  =  Y,hf^S{t-lT),  (9.46) 

/=0 

where  L  is  the  number  of  multipath  components  in  the  channel  response,  T  is  the 
symbol  duration,  and  is  the  complex-valued  channel  coefficient  for  the  /th  path. 
The  channel  coefficients  are  known  at  the  transmitter  and  are  realizations 

of  i.i.d.  zero-mean,  circularly  symmetric  complex  Gaussian  random  variables  with 


variance 


9  Interference  Rejection  and  Management 


241 


E  1^ j  =  2  v/t,  i  and  /.  (9.47) 

It  is  convenient  to  arrange  these  channel  coefficients  for  the  /th  path  in  a  ^  Nt 
matrix  where  i  =  1,2,...,  Nt,  k  =  1,2,  . . . ,  K. 


Fig.  9.16  Tomlinson-Harashima  precoding  applied  to  MIMO  system 


The  MAI  cancellation  is  facilitated  by  use  of  the  QR  decomposition  of  the  chan¬ 
nel  matrix  Thus,  we  express  as 

[H(0)]^  =  QR,  (9.48) 

where  Q  is  an  x  ^  matrix,  such  that  QQ^  =  I,  and  R  is  a  ^  x  ^  upper  triangular 
matrix  with  diagonal  elements  {r//}.  Based  on  this  decomposition  of  the 

signal  to  be  transmitted  is  precoded  with  the  matrix  transformation 

W  =  QA,  (9.49) 

where  A  is  a  ^  x  ^  diagonal  matrix  with  diagonal  elements  l/ra  ,  i  =  1,2,  . . . ,  K . 
The  {r//}  are  real  and  positive  [48].  The  matrix  P  =  pi  is  a  diagonal  K  x  K  matrix 
that  is  used  simply  for  scaling  the  power  of  the  transmitted  signal  and  results  in 
equal  SNR  for  all  users.  Therefore,  we  have  an  effective  channel  matrix  of  the  form 

H(0)wP  =  [QRJ^QAP 

=  pR^A.  (9.50) 

We  note  that  A  is  a  A"  x  K  lower  triangular  matrix  with  unit  diagonal  elements. 
As  a  result,  user  k  sees  multiple-access  interference  from  users  1,2,...,^  —  1.  We 
also  note  that  the  effective  channel  matrix  =  R^A  will  have  full  rank  K, 

provided  that  Nt  >  K. 

By  reducing  this  channel  matrix  to  a  lower  triangular  matrix,  we  can  now  sub¬ 
tract  the  interference  at  the  transmitter  that  each  user  would  normally  observe  at 
their  respective  receivers.  Thus,  when  the  channel  adds  the  same  interference  to  the 
transmitted  signal,  the  received  signal  at  each  receiver  would  be  free  of  interference. 
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Taking  advantage  of  the  lower  triangular  matrix  structure,  successive  interference 
cancellation  is  performed  with  the  feedback  filter  defined  by  the  matrix 

B  =  [I  -  . . . ,  (9.51) 

where  the  matrix  (I— W)  is  used  to  cancel  the  interference  due  to  the  other  users 
that  arises  in  the  current  symbol  interval,  and  the  terms  — , 
are  used  to  cancel  the  interference  due  to  previous  symbols. 

To  ensure  that  the  subtraction  of  the  interference  terms  do  not  result  in  an  increase 
of  transmitter  power,  we  use  the  modulo  operator,  as  in  Tomlinson-Harashima  pre¬ 
coding,  to  limit  the  range  of  the  signal  to  the  boundaries  of  the  signal  constellation. 
Thus,  the  output  of  the  modulo  operators  for  the  ^th  symbol  vector,  as  shown  in 
Fig.  9.16  (for  square  QAM  constellations)  is 

+  Bx^] 

=  -h  Bx„  —  2\/Mz„,  (9.52) 

where  the  modulo  operation  is  performed  on  each  real  and  imaginary  components 
of  the  vector  [s„  +  Bx„],  x^  is  the  ^  x  1  vector  at  the  output  of  the  modulo  operator, 
Sn  is  the  ^  X  1  data  vector,  x„  is  defined  as 

X„  =  [x[  ,  xj_(^_i)]^,  (9.53) 

and  is  a  ^  X 1  vector  with  complex- valued  components  that  take  on  integer  values, 
determined  by  the  constraint  that  the  real  and  imaginary  components  of  x^  fall  in 
the  range  of  \fM).  Therefore,  the  transmitted  signal  vector  is  expressed  as 

s;  =  wpx^ 

=  (9.54) 

and  the  received  signal  vector  is 

L-l 

r„  =  py]H«Wx„_, ■  +  »/„,  (9.55) 

and,  hence, 

P-Ir„  =  x„  +  (H®W  -  I)x„ 

L-l 

+  y]H«Wx„_, •+»?;.  (9.56) 

i  =  \ 

By  substituting  for  B  and  x„  in  (9.56),  it  follows  that 

P“^r„  =  s„  +  1/^  -  2VMz„. 


(9.57) 
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Consequently,  the  MAI  and  ISI  are  cancelled  perfectly,  resulting  in  the  test  statistics 
for  the  ^th  symbol  vector  as  follows: 


y„  = 


(9.58) 


9.4.2. 1  Optimum  Ordering  of  the  Decentralized  Receivers 

The  ordering  of  the  K  decentralized  receivers  affects  the  construction  of  the  K  xNt 
channel  matrix  There  are  K\  possible  column  permutations  of  and, 

hence,  there  is  one  QR  decomposition  associated  with  each  permutation.  In  turn, 
there  are  Kl  transformation  matrices  W  =  QA,  each  of  which  requires  a  different 
transmit  power.  In  order  to  minimize  the  total  transmit  power,  it  is  necessary  to 
search  over  all  the  column  permutations.  Such  an  exhaustive  search  procedure  is 
computationally  time-consuming,  except  for  a  small  number  of  users.  A  method  for 
simplifying  the  search  for  the  optimum  ordering  is  described  in  [19]. 

The  error  rate  performance  of  the  QR  decomposition  method  described  above 
has  been  evaluated  in  [1,2].  Figure  9.17  illustrates  the  symbol  error  probability  as  a 
function  of  the  SNR  (total  transmitted  signal  power  over  all  antennas  divided  by  Nq) 
for  QPSK  modulation,  L  =  1,  2  and  Nt  =  K  =  2.  The  Monte  Carlo  simulation 
results  are  also  illustrated.  The  simulation  results  are  obtained  by  transmitting  1,000 
data  symbols  over  each  of  10,000  channel  realizations. 


Fig.  9.17  Performance  of  optimal-ordered  QR  decomposition  with  Nt  =  K  =  2  and  L  =  1,2 
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Figure  9.18  shows  analytical  results  of  the  symbol  error  rate  performance  for 
QPSK  with  L  =  1  (flat  fading),  K  =  2,  and  Nt  =  2,  3,  4.  We  observe  that  the 
system  performance  improves  with  an  increase  in  the  number  of  transmit  antennas, 
which  reflects  the  benefit  of  spatial  diversity. 


2,3,4 

Figure  9.19  shows  simulation  results  comparing  the  error  rate  performance  of  the 
linear  zero-forcing  (ZF)  and  minimum  MSE  (MMSE)  precoding  methods  with  the 
QR  decomposition  method  for  QPSK  modulation  with  L  =  1  and  K  =  Nt  = 
We  observe  that  the  performance  of  the  QR  decomposition  method  is  better  than 
that  of  the  linear  precoders  at  high  SNR,  but  poorer  at  low  SNR.  It  should  be  noted 
that  the  improvement  in  performance  of  the  QR  decomposition  method  at  high  SNR 
is  obtained  at  a  significantly  higher  computational  complexity  compared  with  the 
linear  precoders. 


9.4.3  Vector  Precoding 

Hochwald  et  al.  [26]  have  described  and  evaluated  the  performance  of  a  vector  pre¬ 
coding  method  in  which  the  data  vector  s  to  be  transmitted  to  the  K  users  is  modified 
by  the  addition  of  a  precoding  vector  with  integer  elements  as  illustrated  in  Fig.  9.20. 
That  is,  the  modified  signal  vector  is  defined  as 

=  s  +  rp. 


(9.59) 
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where  the  components  of  the  signal  vector  are  taken  from  a  square  QAM  sig¬ 
nal  constellation,  r  is  a  real  positive  number,  and  p  is  a  /^-dimensional  vector 
with  complex- valued  elements  whose  real  and  imaginary  components  are  integers. 
Hence,  for  Nt  >  K,  the  transmitted  signal  vector  is 

X  =  Ar  (s  +  rp) ,  (9.60) 


where  the  matrix  A  7  is  given  as 


At 


yU^  (HH^  + 


zero-forcing  criterion 
MMSE  criterion 


(9.61) 


Fig.  9.20  Model  of  MIMO  broadcast  system  employing  vector  precoding 
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The  perturbation  vector  p  is  chosen  to  minimize  the  power  in  the  transmitted 
signal,  i.e., 

p  =  argmin  ||Ar  (s  +  rp)|^  ,  (9.62) 

p 

subject  to  the  constraint  that  the  real  and  imaginary  components  of  each  element  of 
p  are  integers.  Methods  for  solving  the  least-squares  ^-dimensional  integer-lattice 
problem  and  the  selection  of  the  positive  number  r  are  described  by  Hochwald  et  al. 
[26]. 

The  received  signal  vector  is 

r  =  HAr  (s  +  rp)  +  rj,  (9.63) 

where  rj  is  the  additive  noise  vector.  Although  the  perturbation  vector  p  is  not  known 
to  the  receivers,  the  mth  user  assumes  that  its  received  signal  has  the  form 

rm  =  Y  (Sm +rpm)  +  (9.64) 

where  r]'^  includes  the  additive  channel  noise  and  the  MAI  from  other  users  due  to 
the  non-zero  scale  factor  Since  each  user  knows  y  and  r,  the  mth  user  performs 
the  modulo  operation  on  the  real  and  imaginary  components  of  to  remove  pm 
and  passes  the  result  to  its  decoder.  It  is  demonstrated  in  Hochwald  et  al.  [26]  that 
this  vector  perturbation  scheme  achieves  a  signal  diversity  order  comparable  to  the 
diversity  achieved  by  the  maximum-likelihood  detector. 


9.4.4  Lattice  Reduction  Method  for  Precoding 


The  lattice  reduction  method  for  precoding,  described  in  this  section,  is  similar  to 
the  Tomlinson-Harashima  precoding  for  channels  with  known  interference  at  the 
transmitter.  As  in  the  previous  sections,  the  components  of  the  signal  vector  to  be 
transmitted  are  taken  from  a  square  QAM  signal  constellation  and,  thus,  are  ele¬ 
ments  of  a  lattice. 

A  lattice  can  be  expressed  in  terms  of  its  generator  matrix  G  whose  rows  de¬ 
note  a  basis  for  the  lattice,  i.e.,  all  lattice  points  can  be  written  as  a  linear  com¬ 
bination  of  the  rows  of  G  with  integer  weighting  coefficients.  Any  lattice  A  can 
have  many  generator  matrices  and  many  bases  for  representation  of  lattice  points. 
In  particular,  if  F  is  a  square  matrix  with  integer  entries  such  that  det(F)  =  ±1, 
then  F“^  exists  and  its  entries  are  all  integers.  Then,  G'  =  GF  is  a  generator  of 
lattice  A.  The  new  generator  matrix  G'  defines  a  new  basis  for  the  lattice  A.  A 
desirable  property  of  the  modified  lattice  basis  is  that  it  be  an  orthogonal,  or  close- 
to-orthogonal  basis  with  the  lowest  basis  vector  norms.  The  process  of  finding  such 
a  basis  for  a  lattice  is  called  lattice  reduction.  Although  lattice  reduction  in  high 
dimensions  is  an  NP-hard  problem,  a  polynomial-time  suboptimal  lattice  reduction 
algorithm  described  by  Lenstra  et  al.  [29],  called  the  LLL  algorithm,  gives  very 
good  results. 


9  Interference  Rejection  and  Management 


247 


The  lattice  reduction  algorithm  developed  by  Lenstra  et  al.  [29]  was  designed  for 
real  lattices.  Gan  and  Mow  [20]  have  generalized  the  LLL  algorithm  to  lattices  in 
complex  dimensions.  Similar  to  real  lattices,  if  G  is  a  generator  matrix  of  a  complex 
lattice  and  G'  =  GF,  where  F  is  a  square  matrix  with  complex  elements  having 
integer  real  and  imaginary  parts,  such  that  det  (F)  =  ±1  or  det  (F)  =  ±7,  then  G' 
is  also  a  basis  for  the  lattice  generated  by  G. 

We  may  apply  lattice  reduction  to  the  columns  of  H  (HH^)  ^  for  the  zero¬ 
forcing  criterion  or  to  the  columns  of  (HH^  +  ^  for  the  MMSE  criterion. 

Thus,  we  obtain  the  matrix  Fzf  or  Fmse-  For  simplicity,  let  us  consider  the  lattice 
reduction  of  the  matrix  H+  =  (HH^)  ^ ,  which  yields  the  matrices  F  and  F“^ . 

The  precoding  matrix  is  given  as 

Ar  =  H+F.  (9.65) 

The  transmitted  signal  vector  is  x  =  A^s^  where  s'  =  [Fs],  and  the  ele¬ 

ments  of  the  ^-dimensional  signal  vector  s  are  selected  from  an  M-ary  QAM  signal 
constellation.  The  received  signal  vector  is 

y  =  Hx  -h  7/ 

=  HA^s'  +  Tfj 
=  Fs'  +  7/ 

=  mod2y^[s]  -h  jy.  (9.66) 

Thus,  each  receiver  recovers  its  information  symbol  without  interference  from 
other  users. 

The  lattice  reduction  technique  has  been  shown  to  have  a  performance  compara¬ 
ble  to  that  achieved  by  the  vector  precoding  method  described  previously.  In  fact,  the 
signal  diversity  order  achieved  by  the  lattice  reduction  technique  is  comparable  to 
the  signal  diversity  order  obtained  by  maximum-likelihood  detection,  but  the  lattice 
reduction  technique  has  a  much  smaller  computational  complexity.  The  interested 
reader  is  referred  to  the  papers  by  Windpas singer  et  al.  [55]  and  Taherzadeh  et  al. 
[46]. 


9.5  Conclusions 

The  goal  of  this  chapter  was  to  present  an  overview  of  techniques  to  combat  in¬ 
terference  in  wireless  communication  systems.  We  did  this  by  considering  typi¬ 
cal  situations  where  interference  suppression  and/or  mitigation  is  necessary  to  en¬ 
sure  satisfactory  system  operation,  and  these  examples  included  both  uplink  and 
downlink  scenarios.  The  key  waveform  designs  were  single-carrier  direct  sequence, 
multicarrier  direct  sequence,  and  OFDM,  and  included  both  SISO  and  MIMO  sys¬ 
tems.  The  fundamental  conclusion  to  be  drawn  from  the  material  presented  in  this 
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chapter  is  that  many  communication  systems  cannot  function  in  the  absence  of 
explicit  interference  suppression/mitigation  procedures,  over  and  above  the  vari¬ 
ous  modulation  and  error  correction  coding  techniques  that  are  employed  in  their 
design. 
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